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Abstract 
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H next-to-minimal supersymmetric standard model. The vacuum structure of the Higgs potential is 
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I. INTRODUCTION 

Supersymmetric extension of the standard model is one of promising candidates for 
physics beyond the weak scale. In particular, the minimal supersymmetric standard model 
(MSSM) is interesting. However, the MSSM has the so-called fi problem [1]. The /i-term 
is the supersymmetric mass term of Higgs fields. In addition, Higgs fields have soft scalar 
mass term due to supersymmetry (SUSY) breaking. To realize the successful electroweak 
symmetry breaking (EWSB), both sizes of the /i and soft scalar mass must be of the same 
order. Why can the two masses with different sources be of the same order ? That is the 
/i problem. Furthermore, the LHC bounds rrih > 115 GeV [21 E] also give a significant 
constraint on the Higgs sector of the MSSM. 

The next-to minimamal supersymmetric standard model (NMSSM) is the simplest ex- 
tenstion of the MSSM by adding a singlet field 5* |1HT2] ( for a review, see [13]) . The /x-term 
is forbidden by the Z3 discrete symmetry, but it is effectively induced though the coupling 
XSH1H2 in the superpotential after the scalar component of S developes its vacuum expec- 
tation value (vev). Furthermore, such a vev is determined by SUSY breaking terms. Thus, 
the vev is related with the size of SUSY breaking, and the /i-problem can be resolved. 

In addition, the Higgs sector of the NMSSM has significantly different aspects from one 
in the MSSM (see e.g. [I3l|ll]). The lightest Higgs mass at the tree-level can be larger than 
one in the MSSM, because the above coupling term XSH1H2 in the superpotential leads to a 
new quartic term of Hi and H2 in the Higgs potential. The larger value of A would increase 
the lightest Higgs mass. Moreover, the Higgs fields Hi and H2 are mixed with the singlet S 
after the symmetry breaking. Such mixing changes the coupling between the Higgs scalars 
and vector bosons. Thus, the Higgs sector of the NMSSM has a rich structure. 

Furthermore, the behavior of the Higgs potential and its vacuum structure are much more 
comphcated in the NMSSM than those in the MSSM. The Higgs potential of the NMSSM 
should include the realistic minimum, where the successful EWSB is realized. In addition to 
the realistic vacuum, the Higgs potential may include other (local) minima, some of which 
do not break the EW symmetry correctly. If such false vacua are deeper than the realistic 
vacuum, the realistic vacuum may not be realized.^ Hence, our parameter space of SUSY 
breaking terms and dimensionless couplings in the NMSSM would be constrained in order 
to avoid such false vacua such that false vacua are less deeper than the realistic vacuum. In 
fact, numerical studies have been done and also analytical studies along certain directions 
have been carried out [251 - I31] . Then, it was shown that significant parameter regions can 
be excluded by requirement to avoid false vacua. For example, the NMSSMTools^ is the 
famous code to analyze several phenomenological aspects in the NMSSM [35VI37] . and it 
includes some of analytical conditions to avoid certain false vacua. 

Our purpose in this paper is to study phenomenological aspets of the Higgs sector in 
the NMSSM such as the physical mass spectrum of Higgs scalars and their mixing with 



^ In addition, there may be wrong vacua, where squarks and sleptons develop their vevs. On such vacua, 

charge and/or color are broken [HI UHl [T5HM] . 
^ see http://www.th.u-psud.fr/NMHDECAY/nmssmtools.html 



taking into account analytical conditions to avoid false vacua. As false vacua, we consider 
directions, which are new and not included, e.g. in the NMSSMTools. Then, we will 
show that important parameter regions could be excluded by such conditions. Obviously, 
the parameter region, where the Higgs masses are tachyonic on the realistic vacuum, is 
excluded. When the doublet Higgs and singlet Higgs scalars mix sizablely, the lightest Higgs 
boson mass may become tachyonic. When there are such tachyonic modes, there would be 
a wrong vacuum deeper than the realistic vacuum. Thus, the parameter region with the 
tachyonic Higgs boson mass corresponds to the region, where false vacua are deeper than the 
realistic vacuum. Furthermore, the parameter region leading to deeper false vacua would be 
near and outside the region with the tachyonic Higgs boson mass on the realistic vacuum. 
Thus, wider regions of the parameter space, in particular the parameter regions with sizable 
mixing between doublet and singlet Higgs scalars would be excluded by requiring to avoid 
false vacuum. 

This paper is organized as follows. In section [TT| we review the Higgs sector of the 
NMSSM, in particular the realistic vacuum and Higgs boson masses. In section HI, we 
study new false vacua, which can be deeper than the realistic vacuum. In section IV, we 
study numerically implications of our constraints. Section |V] is devoted to conclusions and 
discussions. 



II. REALISTIC VACUUM AND THE HIGGS MASSES IN THE NMSSM. 

We start our discussion with briefly reviewing the realistic vacuum and the masses of the 
Higgs bosons in the NMSSM. The NMSSM is defined by adding a gauge singlet chiral super- 
multiplet S and imposing a global Z3 symmetry to the MSSM. Due to the Z3 symmetry, the 
superpotential consists of only terms involving three chiral supermultiplets, thus dimension- 
ful couplings as a supersymmetric Higgsino mass term and the tad pole term are forbidden. 
In the following, fields with a hat (") symbol represent superfields and those without the 
symbol represent the corresponding scalar fields. The superpotential of the Higgs and the 
singlet superfields is given by 

Wffiggs = -A^^i ■ H2 + ^«:^^ (1) 

where A and k are the Yukawa coupling constants of the Higgs fields, and Hi and H2 are 
the down-type and the up-type Higgs supermultiplets defined as 



H,= "l.H,^Zl (2) 





respectively. 

The realistic vacuum which breaks the electroweak (EW) symmetry successfully can be 
found by minimizing the Higgs potential. The Higgs potential is obtained from F-, Z)-terms 
and the soft SUSY breaking terms given by 



V^oft = mlHlHi + ml^HlH^ + mlS^ S - (xA^SH^ ■ H^ - ^-kA^S^ + h.c^ 



(3) 



where m^^ ^^ g and Aa,k are soft masses and trilinear couplings of the scalars, respectively. 
For the EW symmetry to be successfully broken, the neutral Higgs fields develop vevs while 
vevs of the charged Higgs fields are vanishing. Using the gauge transformations, without 
loss of generality, one can take (-^2'') — ^^^ (-^2) ^ ^^- The condition for vanishing {H^) 
is to require that the charged Higgs scalars have positive mass squareds. Then, the potential 
of the neutral Higgs fields is given by, 

V =\^\S\^ {\Hl\^ + \Hl\^) + \Fs\^ + Vd 

+ m|^ li/iY + rnlM? + "^ll^l' - (AAAi7°i7°^ - ^-kA^S^ + /i.e.), (4) 

where Fs and Vd denote the F-term of S and D-term potential, 

F*s = kS' - \HlHl (5) 

VD = \{9i+gi)m?-\m\'y. (6) 

Here, gi and §2 denote the gauge coupling constants of f/(l)y and SU{2)l, respectively. 
The Higgs sector of the NMSSM is characterized by the following parameters, 

A, K, m]j^, m^2, m|. Ax and A^. (7) 

In the following discussions, we assume that all of soft masses, trilinear couplings and 
Yukawa couplings are real for simplicity. Although the vevs of H^ and S can be complex 
in general under this assumption, it was shown in ^38j that such CP violating extrema are 
maxima rather than minima. Thus, it is reasonable to assume that the neutral Higgs fields 
develop real and non- vanishing vevs while the charged ones do not. Then, we denote vevs 
as 

{H',)=v,, {H',)=V2, {S) = s. (8) 

Furthermore, as was discussed in [27j, the Higgs potential, Q, is invariant under the re- 
placements, A, K, s — !■ —A, — K, — s and A, Vi -^ —A, — vi, therefore we can always take 
A and Vi to be positive while k, /i(= As) and Ax, A^ can have both signs. The exsitence of 
the minima of the Higgs potential is classified according to the signs of k, s and Ax and A^ 



1. for positive k, 

(a) when sign[s] = sign^^A] = — sign^AKJ, the minima always exist. 

(b) when sign[s] = — sign[y4A] = — sign^A^], 

the minima exist if |v4fj| > 3Ai;if2|v4A|/(— |sAa| +fi;|s^|) where the denominator is 
positive. 

(c) when sign [s] = sign [Ax] = sign [A^] , 

the minima exist if \Ai^\ < 3XviV2\Ax\/ {\sAx\ + /t|s^|). 

2. for negative k 



(d) when sign [s] = sign [Ax] = sign [A^] , 

the minima exist if \Ak\ < 3XviV2\Ax\/ {\sAx\ — k\s'^\) where the denominator is 
positive. 

The vevs are determined by the stationary conditions or minimizing the potential, Q, with 
respect to the neutral Higgs fields, 

- — K = X^v COS B(s^ + v"^ sin^ B) — Xkvs^ sin/3 H — gf^t;'^ cos /3 cos 2/3 

+ m1[_^v cos (3 — XAxvs sin (3 = 0, (9a) 

Q = X^v sin /3(s^ + f ^ cos^ (3) — Xkvs^ cos /3 — -.g^v"^ sin /3 cos 2/3 

+ m^2t>sin/3 — A^Afs cos/3 = 0, (9b) 

OT 7" -I 

-- = X^sv"^ + 2k^s^ - Xkv'^s sin 2/3 + m|s - -XAxv^ sin 2/3 + kA^^s^ = 0, (9c) 

where v = \/v\ + f |, tan/3 = ^2/^1 and g^ = gj + gi- The vevs of the doublet Higgs fields 
must satisfy f ~ 174 GeV to give the correct masses to the gauge bosons. Without special 
relations among parameters, when two of the Higgs fields are non- vanishing, the other must 
be non-vanishing, too. Therefore, a non-trivial solution of ^ is as follows; either three Higgs 
fields are non- vanishing or one Higgs field is non- vanishing [52] • This fact is originated from 
the trilinear terms, XAxH^H2S, in the soft SUSY breaking terms and the quartic term, 
Xk,H1H2{S*Y, in the F-term potential. This observation justifies our strategy of analyses 
on false minima of the Higgs potential in the next section. 

It is useful to express the soft SUSY breaking masses in terms of other parameters rewrit- 
ing the stationary conditions, ([9]), 

m^ = — /i^ ^"^1 sin^ j3 m| cos 2j3 + fi { — /i + Ax ) tan /3, 

q^ 2 VA / 



g 

9 9 Z,/ 

m 



/i^ T^mi cos^ /3 H — m| cos 2/3 + /i ( — /i + Ax ) cot /3 

Q'^ 2 VA / 



•) 


(10a) 


J 


(10b) 


jA^f^, 


(10c) 



""' g 

2 2X^ 2 2fi;2 2Xk 2 n X^ Axmi 

mg = -ruz rrrH H ^m^ sm 2/3 + — sm 2/3 

r A^ 5f2 5f^ /i 

where m| = l^f^f^ and /x = As. Thus, given m^, we can use the following parameters, 

A, K, Aa, Ak, tan/3 and /x, (11) 

instead of ([7|). Using these parameters, the minimum of the realistic vacuum, which repro- 
duces the observed Z boson mass, can be written as 

_ ,2 "^|sin^2/3 m|cos22/3 5 , . 

Vmm — -A — —^ h V jjjin, l^iZj 

5 

where V^mjn is the potential involving only s = /i/A, 

vL = Y^ + ^^^^-^^ + x^^l/^'' (13) 

5 



with ?7i| given by (10c). In the following section, we study false vacua and compare their 



depths with (12). 



The mass-squared matrices of the Higgs bosons at tree-level are obtained from (111) by 
expanding the Higgs fields around their vevs. The number of degrees of freedom of the Higgs 
bosons is ten, and three of them are absorbed by gauge bosons via the Higgs mechanism. 
The remaining seven physical degrees correspond to three CP-even Higgs bosons, two CP- 
odd Higgs bosons and one charged Higgs boson. The mass-squared matrix of the CP-even 
Higgs bosons is real-symmetric and denoted as M|j- where i, j runs over 1 to 3 for the 
down-type, the up-type and the singlet Higgs scalars. It is given by 



^h,ii = "^z cos^ j3 + fi i — /i -|- Ax ] tan j3, 
sin^ l3 + fi i — /i + Ax ] cot /9, 



MI22 



K^z 



m, 



4k' 



2 A 

/i + -A,,n + 



A2 Axml 



M, 



h,12 



ML, 



<23 



A2 

2 



2 
2y2A 



nir 



9 
2v^A 



fimz cos (3 



9^ /^ 
sin 2/3 — yU 

V2X 



sin 2/3, 



fi + A; 
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2k 



mz [Ax + -—ft 



firuz Sin (j mz 



X 

2k 



sin/3, 
cos/3. 



(14a) 
(14b) 

(14c) 
(14d) 
(14e) 
(14f) 
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The mass-squared matrix of the CP-odd Higgs bosons, M^, is also real-symmetric and given 

by 



Ml,n 



^1,22 



Ml,, 



2fi 



sin 2/3 

9 
V2X 



Ax + jjJ^ 



-mz 



Ax 



4:K\ ^ 3k , 

— I sm2/3- —A^jj, 



2k 



-/^ 



(15a) 
(15b) 
(15c) 
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where we have removed the Nambu-Goldstone mode. 

It can be shown that the mass of the lightest CP-odd Higgs boson vanishes when k goes to 
zero. This is because the Peccei-Quinn symmetry is restored in this limit. Hence k should 
not be much small to avoid tachyonic CP-odd Higgs bosons. Assuming \fi\ > mz to be 
consistent with non-observation of the charged Higgsinos, we can derive intuitive conditions 



to avoid tachyonic Higgs bosons from (14) and (15). Firstly from (14a), (14b) and (14c), it 



can be understood that the mass of the lightest CP-even Higgs tends to become tachyonic 
if the following conditions are satisfied, 

' K 



^ ij^ + ^'^) ^ 



-m|. 



and/or 



A 



A^fi < -m|. 



(16) 



(17) 



Necessary conditions to avoid the tachyonic Higgs bosons are to require the left-hand side 



of (16) and (17) to be positive, i.e. 



1. For positive k, Axfi and A^fi should be positive. 

2. For negative k, A^fi should be negative and XAxfi > — /t/i^ should be satisfied. 



On the other hand, when k is positive (negative), we can expect from (15b) that Af^jj, should 
be negative (positive) to avoid for the lightest CP-odd Higgs boson to be tachyonic. Thus 
there is a tension between the conditions for non-tachyonic modes of the CP-even and the 
CP-odd Higgs bosons. This tension can be avoided if magnitudes of Ax and A^ should be 
tuned so that both conditions are satisifed. One of the choices of Ax and A^ for positive k 
is obtained as |39| 



^ .2 „„J 4k 2 



Axfi>--fi and — -/i <-A,,fi<0, 



X' 



A2 



K 



(18) 



where tachyonic Higgs bosons can be also avoided when k is negative. The condition (18) is 



enough condition stating that parameters not satisfying this condition result in the tachyonic 
masses of the Higgs bosons. This condition is useful to understand the behavior of the 
tachyonic mass region. As we will show in |IV[ the tachyonic mass region appears near the 
regions given by (18). The conditions imply that the effective Higgsino mass, |/i|, should be 
larger than \Ai^\ . Since /i is of order TeV scale not to introduce the little hierarchy problem, 
\Af^\ must be relatively small. 

Furthermore, we can derive other conditions by taking into account the off-diagonal terms. 
The mass of the lightest CP-even Higgs boson tends to become tachyonic when the mixings 
between the doublet Higgs and the singlet bosons are large. The mixings can be naively 
written as 



2A/i- {\Ax + 2Kfi)sm2(3. 
Requiring this mixing to be vanishing, we have the condition [39] 



A> 



2/i 
sin 2/3 



K 

2-,. 



(19) 



(20) 



Thus the tachyonic Higgs boson tends to appear when the left and right hand sides are not 
comparable. As we mentioned, the supersymmetric Higgsino mass should be about TeV 
scale. This implies that Ax must be taken to be relatively large due to the factor 2/ sin 2/3. 
One might consider that small Ax can be taken when k/X is larger than one so that the 



left-hand side of (20) is tuned. In such parameter region, however, tachyonic Higgs bosons 



do not appear even if Ax is larger than the right-handed side of (20). Indeed, as we will 



see in Sec. IV, the tachyonic Higgs boson appears for small values of \k,/X\. This is because 



the element of the mass-squared, (14c), increases and becomes larger than the off-diagonal 



element (14e) and (14f) as \k,/X\ increases. The mixing of the singlet in the lightest Higgs 
boson is suppressed in this region. 



The mass squared of the charged Higgs boson is 

ml. = ml - ^m| + ^ (^a + jf^) , (21) 

where m^ = ^g2v'^ is the mass squared of the W boson. The charged Higgs boson mass 
squared can also be tachyonic when A is large enough. These mass-squared matrices are 
used in numerical calculations to find tachyonic mass regions. 



III. FALSE VACUA ALONG SPECIFIC DIRECTIONS 

In this section, we show that false vacua in the Higgs potential can be found by considering 
specific directions. Hereafter we analyze the Higgs potential involving only the neutral Higgs 
fields. As discussed in the previous section, when two of the Higgs fields develop their vevs, 
the other must develop its vev to satisfy the stationary conditions. Thus, analyses of the 
Higgs potential are constrained to cases of either one or three non-vanishing Higgs fields. 
Obviously, the vacua, where only one of Higgs fields has non-vanishing vev, are unrealistic, 
and the conditions to avoid such vacua have been studied in [251 132] • The realistic vacuum 
given in dol) is included along the direction where all of three Higgs fields develop their vevs. 
The cases with three non-vanishing Higgs vevs also include false vacua on which the EWSB 
does not occur correctly. Analyses with three non- vanishing Higgs fields are so complicated 
in general that it can not be performed analytically. However, analytical study is possible 
for specific directions in field space along which some of the Higgs fields are related or 
vanishing. In general, minima of the scalar potential appear when positive quartic terms 
balance with negative quadratic and trilinear terms. When the quartic terms in the scalar 
potential, Q, vanish, the false minima appear for large values of the Higgs fields and hence 
these become deeper. Thus we restrict our discussions to three possible cases in which three 
Higgs fields are aligned so that D-term and/or Fg-term are vanishing. In [32], false vacua 
with \Hi\ = \H2\ 7^ and S* 7^ 0, which corresponds to F5 = and Vd = 0, was studied. 
In fact, false vacua deeper than the reahstic vacuum are easily found along these directions. 
Such directions should be avoided to stabilize the realistic minimum. In the following, we 
denote the neutral Higgs fields HI2 as Hi^2 for abbreviation. 

A. Fs = Vd = direction 

Fisrt we consider a direction where both of Fs and Vd are vanishing. The vacua along 
this direction were firstly analyzed in [52] and shown that those can be deeper than the 
realistic vacuum. Along this direction, referred as A, the Higgs fields are not independent 
and related as 

Hi = H2, (22a) 

S^ = -H1H2. (22b) 



Here positive values of Hi and H2 are assumed since we are interested in the false vacua 



near the realistic vacua. Then, k must be positive to satisfy Eq. (22b). When we require 



perturvativity up to the GUT scale (~ 10^^ GeV), i.e. all of the gauge and Yukawa couplings 
to be smaller than \/47r below the GUT scale, the minimum appearing along the direction is 
unacceptable. This is because the Yukawa couplings encounter the Landau pole before the 
GUT scale when tan/3 = f2/'yi = 1 [32]. We will show this numerically in the next section. 
The scalar potential along this direction can be written as 



Va = FH2^ - 2AH2' + m'i/2', 



(23) 



where 



and 



2A3 



K 



A = \^ 



ti\Ay 



\^2\A,\ 



A 



"2 2 I 2 I ^* 2 

m = m„^ +m„^ + -rrig, 

K 



'm 



'H2 



ei = sign[AxHiH2S] = sign [A^] sign [^] ^ 
€2 = sign[K A^S^] = sign [A^] sign [S*] . 



(24a) 
(24b) 
(24c) 



(25a) 
(25b) 



The minimum of the potential becomes deeper when the trilinear term. A, is positive. From 



Eqs. (25), the trilinear term can be always taken to be positive using the sign of S and is 



given as 



A = X^ 



M--^A. 



(26) 



By minimizing the potential of (23) with respect to if 2, the value of if 2 at extremal is 
obtained as 



H- 



2cxt 




(27) 



where rh? < ^-^ is required for if2cxt to be real. Then, the minimum of the potential is 



obtained by inserting (27) into (23) as 



:-^2cxt(^-^2e 



m 



(28) 



To realize the correct EWSB, the following necessary condition is required. 



Va. 



mill _ •'^min- 



(29) 



B. F5 7^ and Vd = direction 



Next, we analyze the direction where Vd is vanishing while Fs is no n- vanishing. We refer 
this direction as the direction B. The vevs of the up-type and down-type Higgs scalars must 
always satisfy the relation 



Hi = H2 



(30) 



for Vd to be vanishing. Similar to III A , the Yukawa couplings blow up before the GUT 
scale, and therefore the minima appearing along the direction should be avoided. 
We parametrize the vev of S as 



S = sign[S]/3H2, 
where P is positive by definition. Then the potential is given by 

Vb = FH2^ - 2AH2^ + m'H^^ 



where 



F 
A- 



1 



2eiA|fi:|/3^ 



U2\\Ax\ - ^e^\n\\A^\(5' 



/3, 



m 



Hi 



+ m^, + m|/3 



IP' 



and 



ei = s\g-n[KH2{S*f] = signf/t], 

€2 = sign[y4Ai^25'] = sign [Aa] sign [S*] , 

eg = sigio.[K A^S^] = sign [k] sign [A«] sign [5] , 



(31) 

(32) 

(33a) 
(33b) 
(33c) 



(34a) 
(34b) 
(34c) 



We can expect that the deepest direction will be found along the positive trilinear terms 
with K > and AxA^ < so that F becomes smaller and A becomes larger. From Eqs. (34), 



the trilinear term can be always taken to be positive, and the quartic and the trilinear terms 
are given by 



F 

A 



1 



2\kI3^ 



\A, 



-kA^P' 



p. 



(35a) 
(35b) 



The extremal value of H2 and the minimum of the potential, VB^mm(/3), along t his d ire ction 
are given in the same form as (27) and (28) by replacing F, A and m? with (35a), (35b) 



and (|33c|). Note that the value of VB,min(/3) depends on /3. Then, the following condition is 

(36) 



required to stabilize the realistic minimum 

VbmM > V, 



10 



for any value of /3. 

It is important to note here that the case with one non- vanishing vev S, i.e., 

Hi=H2 = 0, Sj^O, 



(37) 



is included in the direction B. When we take the limit /3 — )■ oo and -?/2 ~^ by keeping 5* 
to be finite. If the minimum along this direction is the deepest among the ones included in 
the direction B, it can be found automatically by analyzing the direction B. 



C. Fs = and Vd / direction 

The last direction we analyze is that Fs = and Vu ^ 0, refered by the direction C . 

(38) 



From Fs = 0, we have a relation (22b) and parametrize the vev's as 



H^ = aH2, S = sign[S]PH2, 



where a and /3 should satisfy 



/32 = _a, U > 0). 



The potential is given by (23) with 



K 8 



2N2 



A = X 



Ax A^ 

3 



a 



X 



3/2 



- 2 2 2 I 2 I * 2 

m = a rrijj^ + rrijj,^ -\ — amg. 



(39) 

(40a) 
(40b) 
(40c) 



The extremal value of if 2 and the minimum of the potenti al, V c'. min('^) , are given in the 
same form as n27« and (28) by replacing F, A and rh? with (40a), (40b) and (40c). Note 



that the value of Vc',min(«) depends on a. If the extremal values satisfy 
H,l, + H2I, = (1 + a')H2l, = v'^ (174 GeV)^ 



(41) 



this minimum can become the true vacuum. Otherwise, it is a false vacuum. Then, the 
following condition is required to stabilize the realistic minimum 

Vc,min(a) > V"mi„, (42) 

for any value of a. 

The direction C includes the direction along which only H2 has the non-vanishing vev, 
i.e., 



Hi = S = 0, H2^0. 



(43) 
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FIG. 1: Region excluded by the occurrence of a Landau pole on the A-k plane. Solid (red), dashed 
(green), dotted (blue) and dashed-dotted (violet) curves correspond to tan/3 = 5, 3, 1.7 and 1.6, 
respectively. The cut-off of the model is taken as the GUT scale (inside) and 10 TeV (outside). 
The region outside each curve is excluded. 



This can be easily seen by taking a = 0.'^ Furthermore the direction C includes the direction 
A which corresponds to a = 1. Therefore all of the false minima that can be found along 
the directions proposed so far are included in the direction B and C. Hence it is enough to 
analyze the potential along these two directions. 

IV. NUMERICAL ANALYSIS 



In this section, we present numerical results of the constraints from which the false vacua 
studied in the previous section should not be deeper than the realistic one. In addition, we 
take into account that (i) physical masses of the CP-even, odd and charged Higgs scalars are 
non-tachyonic, and (ii) the parameters A, k, and the top Yukawa coupling have no Landau 
pole until the GUT scale (~ 1.6 x 10^^ GeV). For the condition (ii), we solve renormalization 
group (RG) equations at one-loop order [TOl HO] from the EWSB scale to the GUT scale, 
requiring that A and \k\ are smaller than 27r at the GUT scale [131 HI]. 

Figure [T] shows the region excluded by the occurrence of a Landau pole on the X-k, 
plane. We take a cut-off of the model to be the GUT scale (inside) and 10 TeV (outside) 
for references and use the running top quark mass, rrit = 165 GeV as the input. Solid 
(red), dashed (green), dotted (blue) and dashed-dotted (violet) curves correspond to tan/3 = 
5, 3, 1.7 and 1.6, respectively. The region outside each curve is excluded. One can see that, 
for the GUT scale cut-off, A is more constrained as tan/3 is smaller while the upper bound 
on n stays constant around 0.63. This is because RG evolution of A is directly connected 
with the top Yukawa coupling. When tan /3 is small, the top Yukawa coupling at low energy 



The direction with only Hi = can be found by taking a — > and H2 — > 00 keeping Hi finite. 
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No. 


tan^ 


// (GeV) 


Ax (GeV) 


A^ (GeV) 


1 


3 


200 


300 


-50 


2 


3 


200 


-300 


-50 


3 


3 


400 


-300 


-50 


4 


3 


200 


660 


-50 



TABLE I: Parameter sets used in the numerical calculation. 

is large and it grows quickly as the energy scale goes up. Then, A is driven to a large value 
as the top Yukawa coupling grows. On the other hand, RG evolution of k is proportional 
to K^ and depends on the top Yukawa coupling only through A. Therefore, k starts to grow 
after the top Yukawa coupling and A become sufficiently larger than 2tt. As we can see in 
figure [T| the maximum value of A becomes small drastically for tan /3 < 2 and it disappears 
when tan/3 < 1.5. For the 10 TeV cut-off, the upper bounds on A and k do not change with 
respect to tan/3, and are about 1.9 and 1.6, respectively. The result is understood by the 
fact that the RG evolutions of A and n are determined by the values of these couplings at 
the EWSB scale, and the evolutions are almost independent of the top Yukawa coupling. 
This is simply because the cut-off is close to the EWSB scale and the top Yukawa coupling 
does not grow very much even for small tan/3. In the following, we choose moderate values 
of tan /3 to analyze the constraints from the unrealistic minima. 

We use the parameter sets given in Table |T] as illustrating examples. For radiative cor- 
rections to the Higgs potential, we fix the soft masses of the third generation of left-handed 
squark, m^, and right-handed stop and sbottom, m? and rn.? as 



m% = 1000 GeV, mf 



ml 

b 



500 GeV, 



(44) 



respectively. The trilinear term of stop is chosen as nearly maximal mixing so that the 
lightest Higgs boson mass becomes the largest. The point 1 and 4 correspond to the case 
l.(a) and the point 2 and 3 correspond to the case l.(b) explained in Sec|ll| respectively. 
For the point 4, Ax is chosen so that the mixing among the doublet and the singlet Higgs 
vanishes for small k/A. Note that on these points, the minimum is found for positive k. 

In figures [2| we show the contour plot of the lightest CP-even Higgs mass in X-k, plane. 
The values of the Higgs boson mass are shown near each curve or left bottom in the figure. 
The (light blue) filled region represents the tachyonic mass region of all physical Higgs boson 
and the red solid curve represents the Landau pole condition. The figures (a), (b), (c) and 
(d) correspond to the point 1, 2, 3 and 4, respectively. One can see that the tachyonic 
region appears only in small k/A region in figure (a) while it also appears in large k/X 
region in figures (b) and (c). The behavior of the tachyonic mass region can be understood 



intuitively by (18). Firstly, as we discussed in the SecJIIl the mass of the CP-even Higgs 



boson tends to become tachyonic if the diagonal elements of the squared-mass matrix is 



negative. The negative diagonal elements can be avoided when (18) is satisfied for positive 



Ax- The second condition of (18) imposes the lower bound on k/X > \Af,\/An for negative 
Af^. The lower bound is 0.06 or 0.03 for the point 1, 2 or 3, respectively. Hence the diagonal 
elements are negative only for very small k/A. Secondly, the Higgs boson mass also tends 
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FIG. 2: Contour plot of the Higgs masses on the \-k plane. The values of the Higgs masses are 
indicated near each curves and left bottom in GeV unit in the figure. The filled region is excluded 
by the tachyonic Higgs masses and the red solid line represents the Landau pole condition. 



to become tachyonic when the off-diagonal elements are comparable to the diagonal ones. 
The off-diagonal elements become larger as kJ\ becomes smaller. Thus the tachyonic mass 
region appears in the small kJX region for positive A\. For negative A\^ the first condition 



of (18) gives a stronger costraint on k/A because the /i parameter is larger than A^. The 
lower bound on k/A is obtained as |Aa//^| = 1-5 and 0.75 for the point 2 and 3, respectively. 
Thus, the tachyonic mass region appears for relatively large k/A. For the point (d) where 
the parameters are tuned so that the mixing of the doublet and the singlet Higgs bosons 



vanishes, according to (20). In this case, the tachyonic region disappears. This is because 



for small k/A, the second term of (20) is negligible compared to the first term and hence the 
singlet does not mix. The physical mass of the CP-even Higgs bosons are mainly determined 



by the diagonal elements, and the second condition of (18) can not be applied in this case. 



To see the above explanation more concretely, we show the masses of the Higgs bosons 
and the mixings in the lightest Higgs boson for the point 1 in Figures [3j The mixing, 

A^j [i = 1, 2, S), is defined as 



h = N^Hi + N2H2 + NsS, 



(45) 



where hi represents the lightest Higgs boson and J2 ^i — 1- Figures [3J (a) and[3J(b) show 
that the masses of the lightest Higgs boson hi and the second lightest Higgs boson /i2 for 
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FIG. 3: The masses of the hghtest and the second hghtest Higgs boson (top) and the mixings of 
Hi, H2 and S in the hghtest Higgs boson (bottom) for the point 1. The left and right figures 
correspond to k = 0.1 and 0.5, respectively. 



K = 0.1 and 0.5, respectively. Figures [3j(c). and[3j(d) show the mixings of the down-type, 
up-type Higgs Hi, H^ and the singlet S in the lightest Higgs boson corresponding to (a) 
and (b). It is seen in Fig|3}(a) that the mass of hi is the most degenerate to that of /i2 at 
A = 0.2 and it becomes tachyonic at A = 0.4. From FigJ3](c), one can see that for A < 0.2 
the lightest Higgs boson consists of mainly H2 and the mixings of the Higgs bosons are 
constant. The mixings of Hi and S increase at A = 0.2 where the mass of hi is the closest 
to that of /i2- The mixings of Hi and S become maximum at A = 0.4 where the mass of hi 
becomes tachyonic. On the other hand, as is seen in Fig|3](b), the mass of the hi increases 
slowly as A increases. The mixing of Hi and S is small and the main component of hi is 
the up-type Higgs. Thus, the tachyonic mass of hi appears when the mixing of Hi and S 
becomes sizable. 

Figures |4] show the square of the mixing of the up-type Higgs (a) and the singlet (b) in 
the lightest CP-even Higgs boson on A-k plane. From top to bottom, the figures correspond 
to the point from 1 to 3. The values of the mixing are indicated near each curve and the 
filled region is excluded by the tachyonic Higgses. The red curve represents the Landau 
pole condition. From the Figures 4.1. (a) and 4.1.(b), one can see that for the point 1 the 
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FIG. 4: Contour plot of the mixings in the hghtest CP-even Higgs boson on the X-k, plane. 



tachyonic region appears when the squared mixing of H2 (S) in the hghtest Higgs boson 
is smaller (larger) than about 0.5. The qualitative behaviour can be understood by the 
discussion in Sec. |Tlj On the point 1, tachyonic region are mainly determined by the CP- 
even Higgs boson and tachyonic masses of the CP-even Higgs boson appear when k/X is 



small so that (18) and (20) are not satisfied. For small k/X, M^gg becomes very small 



and/or because the mixing with the singlet, MI23! becomes large in this region. The large 
mixing between the up-type Higgs and the singlet results in large mass splitting between 
the lightest and the second lightest Higgs boson, and hence the mass of the lightest Higgs 
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FIG. 5: Excluded region of the lightest CP-even Higgs mass (a) and the coupling of the Z boson 
(b) for the point 1 and 2. 



boson becomes tachyonic. For the points 2 and 3, we can see from Figures 4.2 and 4.3 
that the tachyonic region appears when the squared mixing of the up-type Higgs in the 
CP-even Higgs boson is smaller than 0.7 while the mixing of the singlet is almost zero. 
On these points, the tachyonic mass regions are mainly determined by the CP-odd Higgs 
boson. The diagonal elements of the CP-odd Higgs boson decreases for the negative Ax while 
the off-diagonal element increases negatively. Therefore the physical mass of the CP-odd 
Higgs boson becomes tachyonic unless k/X is relatively large. In this situation, the diagonal 
element M^^^^ decreases and the off-diagonal element M|^2 increases for small k/X region 
due to tan/3 enhancement. Hence the mixing between the up-type and the down-type Higgs 
bosons significantly reduces the lightest Higgs boson mass. The mixing with the singlet does 
not play important role for this choice of the signs because M| 33 is always larger than M| ^^ 
even for negative values. Our numerical results show that the up-type and the down-type 
Higgs doublet mix almost maximally in the lightest Higgs boson, the lightest Higgs boson 
becomes tachyonic. 

Figures |5] and |6] show the constraints from the false vacua on the lightest Higgs mass 
boson (a) and the squared coupling with Z boson normalized by that in the SM (b) in X-k 
plane. Figures [5|l and[5J2 correspond to the point 1 and 2, and Fig|6jl and Figj6|2 to the 
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FIG. 6: The same plots for the point 3 and 4 as the figure [5] 



point 3 and 4, respectively. As the same as in Fig. |2| the values of the lightest Higgs boson 
mass are indicated near the curves in Figures (a) and those of the squared coupling with Z 
boson are indicated in figures (b). The regions excluded by the constraints from the direction 
A, B and C are indicated by filled region coloured as green, orange and purple, respectively. 
The blue filled region represents tachyonic region. In figure [511(a), one can see that the 
constraint from the direction C excludes the region outside the tachyonic region while the 
constraints from the direction B excludes inside the tachyonic region. The region excluded 
by the direction A is very narrow and appears near the horizontal axis. The region excluded 
by the direction B is simply connected from the realistic vaccum along the tachyonic region. 
However, the region excluded by the direction C appears outside the tachyonic region, and 
hence is not connected to the realistic minimum. The region can not be found by only taking 
tachyonic masses into account. From the figure [5J 1(a), one can see that it excludes the large 
region of the lightest Higgs boson mass. From the pj 1(b), one can see that the constraint 
from the C direction excludes the region of the squared coupling from 0.7 to 0.9. This 
implies that the lightest Higgs boson has the SM-like coupling with Z boson and consists of 
almost purely the up-type Higgs in most of allowed region. It is noted that there exists small 
region allowed between the tachyonic mass region and the region excluded by the direction 
C. In this region, the coupling with Z boson is about half of that in the SM. However, 
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FIG. 7: The same figure as figJ5|l(a) and (b). The region of tachyonic squark masses is indicated 
by yellow region. 



as shown in figure [7| squard masses become tachyonic in this region. Thus, this region is 
not reaUstic because colour and charge symmetry are spontaneously broken. In the figures 
5|2(a) and|6}l(a), the constraints from the direction B and C exclude regions outside the 
tachyonic region. The constraint from the direction A also excludes the similar but smaller 
regions. From these figures, large regions on \-k plane are excluded. The allowed region 
corresponds to small A region where the mass of the lightest Higgs boson is smaller than 120 
GeV and the squared coupling is the SM-like. For the point 4, we can see from the figure 
6|2(a) and|6}2(b) that there is no region excluded by the false vacua. As explained before, 



the mixing of the Higgs bosons is vanishing due to the choice of parameters satisfying (20) 



In this case, the false vacua along the direction A, B and C are not deeper than the realistic 
vacuum, and hence the realistic minimum is stable. The lightest Higgs boson mass as well 
as the squared coupling with Z boson are not constrained at all. However, an unnatural 
tuning between A\^ fi and tan/3 is required. The behavior of the constraints or the depth 
of the false minima on the parameters is very complicated and it is difficult to understand 
qualitatively. However, as we have seen, our constraint can exclude sizable region on the 
parameter space which can not be found by tachyonic Higgs mass. Thus, it is important to 
include the constraints from the false vacua in phenomenological studies. 

So far, we have shown numerical resuls for the GUT scale cut-off. The NMSSM with 
the 10 TeV cut-off scale is also interesting in the view of the heavy Higgs boson without 
little hierarchy ^U\ I32H11]- We show the constraints for the point 1 with 10 TeV cut-off 
in figure |8j One can see that the large region on X-k, plane is excluded by our constraints 
although the lightest Higgs boson with the mass 125 GeV and the SM-like coupling can be 
obtained. In the region between the tachyonic Higgs boson mass and the excluded region 
by the C direction, squark masses become tachyonic and hence the region is not allowed 
phenomenologically. 
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FIG. 8: The same figure as figure [SJ 1(a) and 1(b) for the cut-off 10 TeV. 



V. CONCLUSION 

We have studied the mass, the mixings and the couphngs with Z boson of the hghtest 
Higgs boson in the viewpoint of the structure of vacua in the Next-to-Minimal Supersym- 
metric Standard Model. 

In section |TT| we have shown the intuitive conditions for which the tachyonic masses for 
the CP-even and odd Higgs bosons can be avoided. The conditions are derived by requiring 
that the diagonal elements should not be negative and the off-diagonal elements are not 
comparable to the diagonal ones. These give rough bounds on the parameters which are 



useful to understand the behavior of the tachyonic masses. In section III, we have shown 
the new directions, along which unrealistic vacua can appear. The EW symmetry is not 
broken successfully on these unrealistic vacua and hence these vacua should not be chosen 
as our vacuum. We have seen that the depth of the false vacua is characterized by the SUSY 
breaking scale. Therefore the false vacua can become deeper than the realistic vacuum. We 
have also shown that the false vacua studied in the previous works [32] are included in the 
new directions. 



In section IV , we have shown our numerical results on the mass. First we have seen that 



the region of tachyonic Higgs mass appears for small k/A for positive A\ and negative A^ 
while it also appears for large k/A for negative A\ and A^. In the former case, the lightest 
CP-even Higgs boson becomes tachyonic when the mixing of the singlet becomes larger than 
0.5, while in the latter case, tachyonic mass of the CP-odd Higgs boson appears when the 
mixing of the up-type Higgs becomes larger than 0.7. Then, we have shown that the new 
false vacua appear outside the tachyonic Higgs mass region and large k/A region. This 
result implies that the mass and the mixing of the lightest Higgs boson can not be large 
because large value of A is excluded in both cases. In fact, we have shown that by imposing 
the constraint that the realistic vacuum is deeper than the new false vacua, important 
parameter regions for the Higgs mass around 125 GeV can be excluded. The large mixing 
of the up-tye, down-type and the singlet Higgs in the lightest Higgs boson is also excluded 
by the new false vacua, and the lightest Higgs boson consists of mainly the up-type Higgs 
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boson in the allowed region. Then, we have seen that the squared coupling with Z boson of 
the lightest Higgs boson is very close to that of the SM. On the other hand, we have seen 
that the mass and the mixing do not constrain by the tachyonic Higgs mass and the false 



vacua if the parameters satisfy (20). In this case, however, the parameters should be tuned 
so that the mixing of the doublet and the singlet Higgs boson vanishes. 

In general, our constraints exclude wider parameter region than the constraint to avoid 
the tachyonic Higgs masses. In most of cases, the region with small k/A is excluded. Fur- 
thermore, the component of the up-type Higgs in the lightest Higgs boson is also constrained. 
In most of cases, the region with such component less than 0.9 is disfavored. That implies 



that the lightest Higgs boson is SM-like. However, the region with (20) is exceptional, the 
tachyonic modes do not appear and the false vacua are less deep than the realistic vacuum. 
In the end of the conclusion, we comment on the lifetime of the false vacua. As was 
discussed in [32j, if the lifetime of vacua is longer than the age of the universe, a realistic 
vacuum becomes metastable and the parameter space is not constrained. A Euclidean action 
for bounce solutions [151450] should be larger than 400 for the lifetime to be longer than the 
age of the universe. We estimated the Euclidian action for the new false vacuum and 
found the order of 10 - 100. Thus, our results will be still valid if the lifetime is taken into 
consideration. However, detailed studies of the lifetime are important to obtain more serious 
constraints. We will study these aspects in our future works. At any rate, the NMSSM has 
several interesting aspects. It is important to study those aspects of the NMSSM by taking 
into account our new false vacua. 

Acknowledgments 

The authors would like to thank Tevon You for fruitful discussions. T. K. is supported 
in part by a Grant-in-Aid for Scientific Research No. 20540266 and the Grant-in-Aid for 
the Global COE Program "The Next Generation of Physics, Spun from Universality and 
Emergence" from the Ministry of Education, Culture, Sports, Science and Technology of 
Japan. T. S. is a Yukawa Fellow and his work is partially supported by the Yukawa Memo- 
rial Foundation, a Graint-in-Aid for Young Scientists (B) No. 23740190 and the Sasakawa 
Scientific Research Grant from The Japan Science Society. 



[1] J. E. Kim and H. P. Miles, Phys. Lett. B138, 150 (1984). 

[2] 1088223 (2012), 1202.1408. 

[3] S. Chatrchyan et al. (CMS Collaboration) (2012), 1202.1488. 

[4] P. Fayet, Nucl. Phys. B90, 104 (1975). 

[5] P. Fayet, Phys. Lett. B64, 159 (1976). 

[6] P. Fayet, Phys. Lett. B69, 489 (1977). 

[7] P. Fayet, Phys. Lett. B84, 416 (1979). 

[8] H. P. Nilles, M. Srednicki, and D. Wyler, Phys. Lett. B120, 346 (1983). 



21 



[9] J. M. Frere, D. R. T. Jones, and S. Raby, Nucl. Phys. B222, 11 (1983). 
[10] J. P. Derendinger and C. A. Savoy, Nucl. Phys. B237, 307 (1984). 
[11] J. R. Ellis, J. F. Gunion, H. E. Haber, L. Roszkowski, and F. Zwirner, Phys. Rev. D39, 844 

(1989). 
[12] M. Drees, Int. J. Mod. Phys. A4, 3635 (1989). 

[13] U. EUwanger, C. Hugonie, and A. M. Teixeira, Phys. Rept. 496, 1 (2010), 0910.1785. 
[14] U. EUwanger (2011), 1108.0157. 

[15] L. Alvarez-Gaume, J. Polchinski, and M. B. Wise, Nucl. Phys. B221, 495 (1983). 
[16] C. Kounnas, A. B. Lahanas, D. V. Nanopoulos, and M. Quiros, Nucl. Phys. B236, 438 (1984). 
[17] M. Claudson, L. J. Hall, and I. Hinchliffe, Nucl. Phys. B228, 501 (1983). 
[18] M. Drees, M. Cluck, and K. Grassie, Phys. Lett. B157, 164 (1985). 
[19] J. F. Gunion, H. E. Haber, and M. Sher, Nucl. Phys. B306, 1 (1988). 
[20] H. Komatsu, Phys. Lett. B215, 323 (1988). 

[21] G. Gamberini, G. Ridolfi, and F. Zwirner, Nucl. Phys. B331, 331 (1990). 
[22] J. A. Casas, A. Lleyda, and C. Munoz, Nucl. Phys. B471, 3 (1996), hep-ph/9507294. 
[23] T. Kobayashi and T. Shimomura, Phys. Rev. D82, 035008 (2010), 1006.0062. 
[24] Y. Kanehata, T. Kobayashi, Y. Konishi, and T. Shimomura, Phys. Rev. D82, 075018 (2010), 

1008.0593. 
[25] U. EUwanger, M. Rausch de Traubenberg, and C. A. Savoy, Nucl. Phys. B492, 21 (1997), 

hep-ph/9611251. 
[26] K. Funakubo and S. Tao, Prog. Theor. Phys. 113, 821 (2005), hep-ph/0409294. 
[27] D. G. Cerdeno, C. Hugonie, D. E. Lopez-Fogliani, C. Munoz, and A. M. Teixeira, JHEP 12, 

048 (2004), hep-ph/0408102. 
[28] M. Maniatis, A. von Manteuffel, and O. Nachtmann, Eur. Phys. J. C49, 1067 (2007), hep- 

ph/0608314. 
[29] D. G. Cerdeno, E. Gabrielh, D. E. Lopez-Foghani, C. Munoz, and A. M. Teixeira, JCAP 

0706, 008 (2007), hep-ph/0701271. 
[30] R. Franceschini and S. Cori, JHEP 05, 084 (2011), 1005.1070. 

[31] D. G. Cerdeno, J.-H. Huh, M. Peiro, and O. Seto, JCAP 1111, 027 (2011), 1108.0978. 
[32] Y. Kanehata, T. Kobayashi, Y. Konishi, O. Seto, and T. Shimomura, Prog. Theor. Phys. 126, 

1051 (2011), 1103.5109. 
[33] E. Bertuzzo and M. Farina, Phys. Rev. D85, 015011 (2012), 1105.5389. 
[34] K. Hamaguchi, K. Nakayama, and N. Yokozaki (2011), 1111.1601. 

[35] U. EUwanger and C. Hugonie, Comput. Phys. Commun. 175, 290 (2006), hep-ph/0508022. 
[36] U. EUwanger and C. Hugonie, Comput. Phys. Commun. 177, 399 (2007), hep-ph/0612134. 
[37] A. Djouadi et al., JHEP 07, 002 (2008), 0801.4321. 
[38] J. C. Romao, Phys. Lett. B173, 309 (1986). 

[39] U. EUwanger and C. Hugonie, Mod. Phys. Lett. A22, 1581 (2007), hep-ph/0612133. 
[40] N. K. Falck, Z. Phys. C30, 247 (1986). 

[41] D. J. Miller, 2, R. Nevzorov, and P. M. Zerwas, Nucl. Phys. B681, 3 (2004), hep-ph/0304049. 
[42] R. Barbieri, L. J. Hah, Y. Nomura, and V. S. Rychkov, Phys. Rev. D75, 035007 (2007), 

hep-ph/0607332. 

22 



[43] P. Lodone, Int.J.Mod.Phys. A26, 4053 (2011), 1105.5248. 

[44] E. Bertuzzo and M. Farina (2011), 1112.2190. 

[45] S. R. Coleman, Phys. Rev. D15, 2929 (1977). 

[46] S. R. Coleman, [Erratum-ibid] Phys. Rev. D16, 1248 (1977). 

[47] C. G. Callan, Jr. and S. R. Coleman, Phys. Rev. D16, 1762 (1977). 

[48] S. R. Coleman and F. De Luccia, Phys. Rev. D21, 3305 (1980). 

[49] K.-M. Lee and E. J. Weinberg, Nucl. Phys. B267, 181 (1986). 

[50] M. J. Duncan and L. G. Jensen, Phys. Lett. B291, 109 (1992). 



23 



